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Abstract. - Molecular motors transduce chemical energy obtained from hydrolizing ATP into 
mechanical work exerted against an external force. We calculate their efficiency at maximum 
power output for two simple generic models and show that the qualitative behaviour depends 
crucially on the position of the transition state or, equivalently, on the load distribution factor. 
Specifically, we find a transition state near the initial state (sometimes characterized as a "power 
stroke") to be most favorable with respect to both high power output and high efficiency at 
maximum power. In this regime, driving the motor further out of equilibrium by applying higher 
chemical potential differences can even, counter-intuitively, increase the efficiency. 
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Introduction. — Molecular motors are essential for 
directed transport within the cell [1]. They typically op- 
erate under nonequilibrium conditions due to the unbal- 
anced chemical potentials of molecules like ATP or ADP 
involved in the chemical reactions accompanying the mo- 
tor steps. In contrast to macroscopic engines, fluctua- 
tion effects arc important thus allowing for backward steps 
even in directed motion. The stochastic dynamics of these 
motors under an applied load force can be probed experi- 
mentally by single molecule assays (see, e. g., for kinesin 
[2], myosin [3-5] or ATPase [6]). Generically, such biomo- 
tors are modelled either in terms of continuous "flashing 
ratchets" [7-10] or by a (chemical) master equation on a 
discrete state space [11-17]. 

For macroscopic engines working between two heat 
, baths at temperatures T2 > Ti , efficiency is bounded by 
' the Carnot limit rjc = 1 — Ti/T2. Since this limit can only 
, be achieved by driving the engine infinitcsimally slowly, 
thus leading to an infinitcsimally small power output, it is 
arguably more meaningful to characterize engines by their 
efficiency at maximum power [18,19]. This quantity has 
been studied for more than 30 years under the label of 
"finite-time thermodynamics" [18-22]. Recently, this con- 
cept has been transferred to microscopic (Brownian) heat 
engines in a variety of different model systems [23-25]. 

In contrast to heat engines, biomotors are driven by 
chemical potential differences. The efficiency of such mo- 
tors is bounded by jymax = 1 [26]. This bound can only 
be reached in an equilibrium situation corresponding to 
a vanishing power output of the motor. In analogy with 



heat engines, we here propose to investigate such motors 
under the condition of maximum power output. 

We start with a simple model system for a chemically 
driven biomotor [11] and show that the qualitative results 
also apply to a more realistic motor model involving a 
second cycle. In this second model, a futile cycle leads 
to ATP consumption even at the stall force. Mechanical 
and chemical cycles arc no longer tightly coupled. In both 
cases, the efficiency at maximum power crucially depends 
on the position of the transition state or, equivalently, on 
the load distribution factor. In fact, a transition state 
near the initial position is most favorable with respect to 
a maximal motor power output. For the efficiency at max- 
imum power, we obtain two counter-intuitive results : (i) 
it increases when the transition state position is changed 
in such a way that the power output rises and (ii) it can 
increase when the system is driven further out of equilib- 
rium by a higher chemical potential difference. 

Model I. — We first consider a linear molecular 
motor with equivalent discrete states (sites) Xn {n = 
0, ±1, ±2, . . . ) with distance / between the sites and next- 
neighbour transitions between these states subject to a 
force F in backward direction, see Fig 1. Forward re- 
actions are assumed to be driven by ATP molecules 
with chemical potential hatp and backward transitions 
by ADP and P molecules with chemical potentials //adp 
and /ip, respectively, 



ATP + Xn^ Xn+1 + ADP + p. 



(1) 
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Fig. 1: (Color online) Scheme of the free energy landscape V{x) 
of a molecular motor with step size I. The transition state 
position S determines the force dependence of the transition 
rates. 



If the dilution of all involved species is high, wc can as- 
sume mass action law kinetics for the rate constants. Ad- 
ditionally, we assume the usual force dependence of rate 
constants [11] such that the transition rates for forward 
and backward steps are given by 



w 



CADPCp/fc-e«i-^)^\ 



respectively. Here, are the (dimensionless) concentra- 
tions of i molecules {i = {ATP, ADP, P} ) and P = 1 //sbT 
with Boltzmann's constant fcs- The bare reaction rates 
fc+, are concentration independent. The load distri- 
bution factor < S < 1 characterizes the location of the 
transition state, see Fig. 1. It can vary between the ex- 
treme cases (5 = (sometimes characterized as a "power 
stroke" [27] ) and S = 1, where forward or backward rate 
constants, respectively, no longer depend on the force. The 
chemical potential of the involved molecules is 



(4) 



with a reference value fi^. Thermodynamic consistency 
requires 



/w 



,(AM-F;)/(fcBT) 



(5) 



where Afj. = /iatp — Madp — yUp. The (mean) velocity of 
the motor can then be calculated as 



-S)Fl 



V = l{w^ — W ) ~ k CADPCp/ 

(6) 

Thermodynamic quantities for each single transition can 
now be defined [15,28,29] on the basis of the transition 
rates. The chemical work applied during one forward step 
is just the chemical potential difference Wchem = A^. The 
mechanical work delivered by the molecular motor during 
a single forward step against the applied force F is given 
hy W = Fl. Since all states are equal, the internal en- 
ergy does not change, AE = 0, and thus the difference 
Q = VFchom — W is dissipated as heat in the thermal envi- 
ronment. The efficiency rj of this chemical motor is given 




Fig. 2: (Color online) Maximum power W* in units of 
fc^CADPCpfceT' exp(/3A/i) as a function of the position of the 
transition state S. 



(2) . . 

by the ratio of mechanical work and chemical work applied 

(3) by the chemical potential difference [26] as 



W 



V 



Fl 



(7) 



With the force velocity relationship (6), the power output 
follows as 



W = Fv = k-CADpcplF 



^P(Ati-SFl) _ g/3(l-5)F/ 



(8) 



The power W is zero for F ^ Q. When the force ap- 
proaches the stall force F F^^ = A/i//, where the ve- 
locity vanishes, the power output also becomes infinites- 
imally small. Thus, there is an optimal force F*, where 
the power output is maximal for a given chemical potential 
difference A/i. This optimal force is given by dW/dF = 
which leads to the implicit relation 



^0Au , , _.01F' 1 + (1 - S)[31F* 
1 - SplF* 



(9) 



For given /3A^, the scaled optimal force (31F* depends 
only on the parameter i5. The optimal power W* is shown 
in Fig. 2 as a function of the transition state position 5. 
Numerical results for the efficiency at maximum power rj* 
are shown in Fig. 3. Note that the latter results are quite 
universal since no kinetic parameters enter these graphs. 
Both, power output and efficiency increase with decreas- 
ing 5 and thus, a transition state near the initial position 
{5 = 0) is most favorable. Previously, it has been specu- 
lated [27] that such a mechanism, where forward rates are 
almost independent of the force, is realized in molecular 
motors in order to reach a large motor velocity (corre- 
sponding to a high power output). Beyond corroborating 
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Fig. 3: (Color online) Efficiency at maximum power r/* as a function of (a) the position of the transition state 5 and (b) the 
chemical potential difference A/x for Model I. 



this idea, we find as a new result that small 5 also leads 
to a higher motor efficiency at maximum power. This 
is somewhat counter-intuitive since an increase in power 
usually leads to a decrease in efficiency. 

In the limit of small chemical potential differences 
(where the motor works in a linear response regime near 
equilibrium), efficiencies at maximum power can be ob- 
tained analytically. In this limit, the stall force also be- 
comes small and thus the exponentials in (6) can be ex- 
panded and truncated after the first order in A/i and F 
leading to the approximate force-velocity relation 



V Ki k CAi)pcpf3l{A^ — Fl). 



(10) 



In analogy to the linear response result for heat engines 
[30] , the efficiency at maximum power universally becomes 
77* = 1/2. Beyond linear response, as a somewhat surpris- 
ing result, the efficiency at maximum power increases for 
increasing chemical potential differences for positions of 
the transition state S < 1/2, compare Fig. 3b. Usually, 
dissipative cost increases when the system is driven further 
out of equilibrium. 

Model II. — In order to check the generality of the 
results obtained for the (simple) Model I, we now calcu- 
late the efficiency at maximum power for a more involved 
motor model. Recent experiments focussing on the back- 
steps of kinesin [2] indicate that a realistic motor model 
should comprise at least one additional cycle [2, 16, 17] 
leading to non-zero dissipation even at stall force. Such a 
mechanism with additional motor cycles presumably also 
applies to myosin motors which have a similar molecular 
structure [31]. In order to capture the main experimental 
finding of ATP-driven backsteps from Ref . [2] , we propose 
a minimal model as shown in Fig. 4. This model is very 



similar to recent kinesin models [2, 16, 17] and thus cap- 
tures experimental findings qualitatively. Binding and hy- 
drolyzing ATP leads to the unbinding of one motor head. 
The elastic energy then leads to a biased diffusive search 
for the next binding site. For high load forces, the prob- 
ability of a backstep increases. Note that such backsteps 
involve ATP consumption and thus decrease the coupling 
ratio between chemical and mechanical motor cycles. The 
force dependence is modelled as 



H\ ~ "'21'^ J "'12 ~ 1^12^ 



21 ~ f^21^^ " I ^12 ~ ^12^^^ 



(11) 



with the transition state located at 61^2 for forward 
and backward steps, respectively. Thermodynamic con- 
sistency requires fc^i/fcj~2 = exp(/3A£'), fc2i/^']*2 ~ 
exp(/3A_E), and 'W12/W21 = exp[/3(A/i — AE')] where IS.E is 
the potential energy difference between state 1 and state 
2. Given all rate constants, the steady state can be calcu- 
lated as 



Pi 

P2 



W2I + + W21 



W2I + W21 + W21 + Wl2 -\- W12 + 1^12 



, (12) 



W12 + W{2 + W12 



W21 + W^i +W21+ W12 + W12 + W12 



(13) 



with motor velocity 

V = I [p{{wt2 - W^2) +P'2{wtl - W'21)] • (14) 

The power output of the motor is M'^ = Fv. Chem- 
ical work is applied to the motor only in the (vertical) 
transitions involving ATP and ADP -\- P. In such a step, 
chemical energy of amount A/i is transferred to the system 
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Fig. 4: (Color online) Two-cycle model for a molecular motor (Model II). (a) Scheme of the reaction pathways, (b) Definition 
of rate constants, (c) Efficiency at maximum power as a function of the transition state position 5 = Si — §2 with ^121 = 
0.275/s, fcai = 0.7/s, k^^ = 1.8/s, AE = fOfcsT. (d) Relation between force [in pN] and velocity [in nm/s] for Si = 0.004, S2 = 
0.024, W21 = 0.275/s, fcJi = 0.7/s, fe+ = 1.8/s, A£ = lOfcijr, = 20kBT,l = 9{kBT)/pN ~ 36nm compared to corresponding 
data from myosin experiments [5]. For these parameters, the optimal force is F* ~ 1.5pA'^ and the efficiency at maximum power 
becomes 77* ~ 0.18. 



and thus the chemical work per unit time is 

PK;hcm = A/.t(p-Jwi2 -_P2^2l)- (15) 

The motor efficiency rj = W/Wchem can then be calculated 
for a given set of rate constants and a given force. 

We again ask for the optimal force leading to a maxi- 
mal power output. We recover the qualitative results of 
Model I for the maximum power (data not shown) and 
the efficiency at maximum power, see Fig. 4c, also in this 
(more realistic) model of a molecular motor. Specifically, 
the largest efficiency can be achieved for both transition 
states near the initial position {S = Si = 0). For 

small 6, the efficiency first increases with increasing chem- 
ical potential difference Afi until it reaches a maximum. 
The advantage of a transition state near the initial posi- 
tion with respect to high power output and high efficiency 
thus seems to be a quite general characteristics for molec- 
ular motors. Note that efficiencies are generally lower due 
to the ATP-driven backward steps leading to additional 
dissipation in such models with additional cycles. 

Discussion. — In summary, we have first investigated 
a simple genuine model of a molecular motor under the 
condition of maximum power output. As our main result, 
we find that a transition state near the initial position 
yields both the largest power output and the largest effi- 
ciency at maximum power. Qualitatively, this behaviour 



is also recovered in a more realistic model involving a sec- 
ond motor cycle. The advantage of a small load distri- 
bution factor S with respect to large power output and 
high efficiency at maximum power thus seems to be quite 
generic. We have assumed < 5 < 1 throughout this 
letter. Recently, it has been proposed to use 5 < to in- 
terpret non-monotonous force velocity relations [32]. For 
negative (5 < 0, the efficiency at maximum power is even 
larger than for (5 = in both our model systems. 

For both models, the efficiency at maximum power can 
increase when the system is driven further out of equilib- 
rium by a higher chemical potential difference. This result 
should be distinguished from previous work predicting a 
maximal efficiency for non-zero chemical potential differ- 
ences [17,26]. In our first model, the efficiency decreases 
monotonically as a function of the chemical potential dif- 
ference for a given load F. 

For kincsin motors, it is difficult to find clear evidence 
for a putative design principle of a transition state near the 
initial position. While previous studies have found small 
(5 < 0.1 for the main motor step [13], a transition state in 
the range 6 ~ 0.3. ..0.65 has been extracted recently [16]. 
For myosin motors, small S has been reported for the main 
motor step {6 < 0.1 [4], [14]) , which is also supported by 
our model II as detailed in the following. 

In Fig. 4d, the force velocity relation of our model with 
appropriate rate constants is compared to a recent myosin 
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experiment [5]. For forces F > 0, the experimental data, 
including a step at ~ ^^J■/l, is captured by our sim- 
ple model. For negative forces, our model shows a second 
step at ~ —Afi/l which is not present in the expcrimctal 
data. Both steps in the theoretical model are explained by 
the fact that for large forces |F| ^ A^//, vertical (ATP 
driven) transitions are slow compared to the horizontal 
(force dependent) transitions and the motor is basically 
pulled by the load. The step at negative forces could be 
eliminated by introducing a force dependent energy dif- 
ference AE{F) which corresponds to an additional force 
dependence of the vertical transitions. However, instead 
of introducing new parameters, we keep the simple model 
II and extract from the crude fit shown in Fig. 4d transi- 
tion state positions i5i = 0.004, 62 = 0.024. These values 
mainly determine the behaviour of the molecular motor at 
large loads |F| ^ A^/l. With these parameters, we find 
an efficiency at maximum power of ry* ~ 0.18 which, given 
the strong decay of ry* with increasing S, is quite close to 
the optimal value rj* ~ 0.31 for 5i = 62 = 0, see Fig. 4c. 

We do not claim that our simple model can explain all 
aspects of myosin motility. Rather, we have chosen model 
II in order to probe the robustness of our main results 
concerning the efficiency at maximum power. In order to 
construct a comprehensive model of myosin, more experi- 
mental data seem to be necessary. 

If future experiments confirm the indication that myosin 
motors have a transition state near the initial position for 
the main motor step (corresponding to an almost force- 
independent forward rate), it would be tempting to specu- 
late whether evolutionary pressure for efficiency and large 
power has led to this characteristic. For a more com- 
prehensive answer to this question, however, other evolu- 
tionary goals like speed, robustness, and high processivity 
should be considered. Likewise, the dependence of our 
results on the interaction between single motor domains 
needs to be explored in future work. 
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